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ABSTRACT 

The Advanced Liquid-Mirror Probe of Asteroids, Cosmology and Astrophysics (AL- 
PACA) is a proposed 8-meter liquid mirror telescope surveying ~ 1000 deg 2 of the 
southern-hemisphere sky. It will be a remarkably simple and inexpensive telescope, 
that nonetheless will deliver a powerful sample of optical data for studying dark en- 
ergy. The bulk of the cosmological data consists of nightly, high signal-to-noise, multi- 
band light curves of SN la. At the end of the three-years run ALPACA is expected to 
collect > 100, 000 SNe la up to z ~ 1. This will allow us to reduce present systematic 
uncertainties affecting the standard-candle relation. The survey will also provide sev- 
eral other datasets such as the detection of baryon acoustic oscillations in the matter 
power spectrum and shear weak lensing measurements. In this preliminary analysis 
we forecast constraints on dark energy parameters from SN la and baryon acoustic os- 
cillations. The combination of these two datasets will provide competitive constraints 
on the dark energy parameters under minimal prior assumptions. Further studies are 
needed to address the accuracy of weak lensing measurements. 
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1 INTRODUCTION 



Over the past decade a picture of the Universe has emerged 
from a number of cosmological observations which have 
shed light o n its geometry, matter c o ntent and clustering 
properties JPe Bernardis et al. 20001: iPercival et al. 200 fl 
ISoergel et al. 20031 : iTeemark et al. 20041) . Most astonishing 
amongst these findings is certainly the existence of "dark en- 
ergy" which dominates the total energy density of the Uni- 
verse and which i s responsible for its present state of acce l- 
erated expansion jRiess et al. 19981: iPerlmutter et al. 19991) . 

Despite the accurate measurements of the Cosmic Mi- 
crowave Background anisotropies, the SN la luminosity dis- 
tance observations and the detailed mapping of the large 
scale matter distribution, the nature of this exotic compo- 
nent still remains unknown and little progress has been made 
so far in this field, both at theoretical and observational 
level. 

The cosmological constant has been advocated as the 
simplest candidate, since it is physically motivated repre- 
senting the energy contribution of the vacuum. However 
we have no convincing explanation as to why the observed 
value is extremely small compared to particle physics expec- 
tations. Alternative scenarios have flourished in the recent 



past, although none seems to provide a consistent particle 
physics formulation. For a recent review see Padmanabhan 
(2005). 

Supernova type la observa tions give the str ongest evi- 
dence in favor of dark energy jRiess et al. 2 004). nonethe- 
less the limited number of measurements are still affected by 
large experimental and systematic uncertainties which pre- 
vent further insights. Consequently many questions about 
the properties of dark energy still remain unanswered. For 
instance it would be crucial to know whether the dark en- 
ergy is time dependent or not, whether it clusters, whether 
its effects are caused by an exotic form of matter or by a dif- 
ferent behavior of gravity on the large scales. Recent data 
analysis have shown that current observations lack the nec- 
essary accuracy to address these very issues (see for instance 
Corasaniti et al. 2004 and reference therein). 

Complementary to the SN la observations, sev- 
eral other tests of dark energy have been stud- 
ied in a vast literature including baryon acous- 
tic oscillations in the galaxy power s pectrum 
jBlake fc Glazebrook 20031: ISeo fc Eisenstein 20031). cluste r 
number counts JWeller fc Battve 20031 IWang et al. 20041) 
weak lensing shear measur ements llHu^^am^OM : 
iBenabed fc Van Waerbeke 20041) just to mention a few. 
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Improving the accuracy and reducing the systematics is 
the necessary condition for these methods to be effective. 
However they are all limited to some extent by degeneracies 
between dark energy and the other cosmological param- 
eters. It is therefore unlikely that one single experiment 
may reveal the nature of dark energy. Progress seems only 
possible through a synergy of these different tests. 

The Advanced Liquid-mirror Probe of Asteroids, Cos- 
mology and Astrophysics (ALPACA) is a proposal for the 
realization of a ground based liquid-mirror telescope sur- 
veying a strip of sky of the southern-hemisphere. The survey 
will collect a huge amount of data divided in several datasets 
pertaining to SN la light curves, high redshift galaxies and 
weak lensing. 

In this paper we study the constraints on dark energy 
from the estimated distribution of SN la observed during 
the nominal duration of the survey, and from the expected 
level of detection of the baryon acoustic oscillations in the 
galaxy power spectrum. The paper is organized as follows. 
In Section[5]we give an overview of the project. In Section^] 
we discuss datasets that the survey is expected to provide, 
in particular SN la observations, baryon acoustic oscillations 
and weak lensing. In Section 2] we introduce the Fisher ma- 
trix analysis and in Section|^|we present our results. Finally 
in Section||j]we discuss our conclusions. 



2 ALPACA PROJECT 

ALPACA is an 8-meter liquid mirror telescope which will in- 
corporate a 3°-wide imaging field and will be sited on Cerro 
Tololo (CTIO) in northern Chile. It is amazingly cost effec- 
tive, delivering more information than virtually any compet- 
ing survey instrument but at a small fraction of the cost. It 
will be a remarkably simple and inexpensive telescope that 
nonetheless will deliver powerful samples of optical data that 
are suitable for dark energy studies. 

The concept of large, liquid-mirror telescope (LMT) is 
fairly new, although LMTs have successfully operated for 
over a decade (Cabanac, Borra & Beauchemin 1998; Hick- 
son & Mulrooney 1998). Serious difficulties which compro- 
mised the past optical performance of LMT have now been 
eliminated. The result of these advancements in the liquid- 
mirror technology is seeing-dominated imaging. The next 
generation of LMT such as ALPACA can therefore compete 
with standard large optical telescopes for many projects. 

ALPACA imaging quality is expected to be limited only 
by the seeing at CTIO, hence < 1 . The cosmological sur- 
vey field will cover 775 deg 2 in five bands (it, b, r, i, z) with 
limiting nightly AB magnitudes in the range 23rd to 25th 
over a large fraction of the field. Four of the five photo- 
metric bands u (310-410nm), b (415-555nm), r (560-745nm) 
and i (750-1050nm) are spaced evenly in log(X) (A being 
the wavelength), so that K-correction uncertainty vanishes 
at the special ratios of (1 + z) — 1.25, 1.58 and 1.98. A 
fifth band z (950-1050nm) is especially useful for isolating 
high-redshift objects. An additional 100 — 200 deg 2 will be 
covered at lower sensitivity. ALPACA will observe the field 
passing directly overhead at CTIO, some 3° across, of which 
2.5° are well suited for cosmological work. The focal plane's 
CCDs operate in drift-scan mode, each CCD reobserves the 
same strip of sky each night, making image subtraction and 



further analysis procedures stable, simple and redundant. 
We are also exploring installation of a multiobject spectro- 
graph. The nominal duration of the survey is 3 years, at 
which time the flux depth will approach the crowding limit. 
The total amount of data collected during the three years 
run is expected to be of several petabytes. 

Observations of SN la standard candle luminosity dis- 
tance are expected to provide the bulk of the cosmological 
data, other datasets will include weak lensing, baryon oscil- 
lations, cluster number counts, strong-lens time delays, and 
cross-correlation multiple-waveband methods, such as ISW- 
correlation with CMB maps and X-ray galaxy clusters. 



3 DATA PRODUCTS 

3.1 SN la standard candles 

Supernovae Type la light curves are related by one pa- 
rameter family relati on between the time-width and the 
maximum-light peak JPhillips 19931 : lHamuv et al. 1996ft . As 
consequence of this SN la are standardizable candles and 
can be used as distance indicators (Riess, Press & Kirsh- 
ner 1996; Perlmutter et al. 1997; Phillips 1999). Models of 
supernova have been intensely investigated for decades, but 
only in the past few years have three-dimensional numerical 
simulations provided a better understanding of the physical 
processes which take place during the explosive phase (see 
Hoflich et al. (2003) for a recent review). 

The most accredited scenario consists of a car- 
bon/oxygen (C/O) White Dwarf (WD) accreting material 
from an evolved companion star. As the WD reaches the 
Chandrasekar mass, the compressional heating triggers a 
thermo-nuclear explosion in the WD's core that propagates 
to the outer layers. However we still lack of a comprensive 
understanding of the fine details of the explosive phase which 
has limited the possibility to make robust quantitative pre- 
dictions for the SNe light curves and spectra. In fact there is 
evidence that a number of secondary effects are responsible 
for deviations of the SN light curves from the one parameter 
family relation. These may depend on conditions prior to the 
explosion, such as the initial metallicity of the WD, the C/O 
ratio as function of the mass of the progenitors, and features 
of the explosive phase, such as rotational support, magnetic 
fields, convection structure in the deflagration front, view- 
ing angle, ejecta asymmetries, effects of the companion on 
the ejecta, circumstellar interactions just to mention a few. 
Some of these may be purely stochastic, while others may 
leave clues in multiband lightcurves which will allow their 
induced variance to be removed or reduced. Several theoret- 
ical treatments have attempted to predict the effect on light 
curves of these various intrinsic factors (Hoflich, Wheeler & 
Thielemann 1998; Dominguez & Hoflich 2000; Mazzali et al. 
2001; Pinto & Eastman 2001; Timmes et al. 2003). Unfor- 
tunately, they tend not to agree well on the size, sign or 
even the rough nature of these effects. In addition extrinsic 
effects, such as varying extinction laws or gravitational lens- 
ing may also influence the luminosity and other lightcurve 
observable as well. 

The use of the brighter-slower relation has already al- 
lowed for a reduction of the dispersion in the magnitude- 
luminosity distance to ~ 15% r.m.s. (Phillips 1993, 1999; 
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Table 1. The redshift distribution of SN la expected from ALPACA, where z corresponds to the center of 
the bin. 
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Riess, Press & Kirshner 1996; Perlmutter et al. 1997; Gold- 
haber et al. 2001). Several studies have shown that this 
scatter can be further reduced using color data. For in- 



stance using the CMAGIC method (IWang et al. 20031) to 
find when B — V color reaches a certain level after maxi- 
mum light and thereby inferring luminosity, the dispersion 
can be limited to 8 — 10%. This method still requires an accu- 
rate measurement of the time-width of the maximum-light 
peak. On the contrary using the B — V color 12 days after 
maxi mum- light (AC12 ) can reduce the scatter to ~ 7% or 
less (IWang et al. 2005T) without the need of the maximum- 
light peak width. 

The ALPACA survey strategy is optimal for detecting 
a huge sample of supernovae, providing superlative qual- 
ity dataset of nightly, high signal-to-noise, multiband SN la 
light curves. At the end of the three years run ALPACA 
is expected to observe > 100, 000 supernovae distributed in 
ten equally spaced redshift bins in the range 0.2 < z < 1, 
and several hundred low-z supernova (z < 0.2). Hence a pri- 
mary utility of the ALPACA dataset would be to improve 
the accuracy of the SN la standard candle relation and re- 
duce systematic effects. In Table we list the expected SN 
la redshift distribution. 1 

As discussed above, the theoretical uncertainties in the 
physics of SN la leave us without indication on how indi- 
vidual characteristics of the progenitors or peculiar condi- 
tions during the explosive phase may systematically affect 
the standard-candle relation. If such effects exist we can 
only hope to identify them through the component anal- 
ysis of large sample of nightly, high S/N, multiband light 
curves. Although it is not a priori obvious which observa- 
tional features might tie to luminosity, we might already 
have hints. For example, the CMAGIC technique (Wang 
et al. 2003) uses intensively sampled multiband lightcurves 
to produce a color-based correction to the maximum-light 
standard-candle relation. However there is still more infor- 
mation in such light curves. As an example, B-band "bump" 
of additional luminosity above the observed, or linear Mb 
versus B — V evolution can add up to ~ 0.5 mag to the B lu- 
minosity near maximum light, which indeed seems to be at 
least roughly correlated with the SN peak luminosity. Such 
features might become manifest as additional higher order 
components to the brighter-slower relation. Using our high 
signal-to-noise, densely sampled, multicolor lightcurves, we 
expect to be able to measure from the component analysis 
even 4th and 5th-moments, and their possible correlations 
with luminosity. 

Furthermore multiband lightcurve shape may be tied to 
the overall SN luminosity. For instance, it has been suggested 
that under multiple scattering conditions, light echoes can 
affect the effective luminosity and color of SNe la events 
(Patat 2005). Single scattering has negligible effect at the 



We thank Ben Johnsonn for computing the expected SN la 
redshift distribution. 



level of accuracy required in the SN la standard-candle rela- 
tion, but in the case of optical-thickness of order of unity, cir- 
cumstellar echoes may prove important even for SNe la near 
maximum light, to the point of affecting stretch or color- 
based luminosity corrections significantly. Such an echo will 
affect various bands differently, and may be or may not be 
associated with high extinction, depending on the patchi- 
ness of the surrounding dust. These characteristics makes 
ALPACA nightly multiband lightcurves an optimal bench- 
mark to measure this effect. 

In addition the large ALPACA dataset will allow us 
to search for possible subclasses of SN la and to the dis- 
covery of new standard-candle relations at different red- 
shifts. Since we can already disentangle SN la from core- 
collapse SN via their multiband light curves alone (Johnson 
and Crotts 2005), we can simply rely on spectroscopy of 
a few representative members for characterization of each 
la subclass we discover (not necessarily for type classifica- 
tion of every SN). In fact hydrogen-rich supernova (Type 
II) can be distinguished from the hydrogen-poor (la, lb and 
Ic) through the UV deficit, since UV radiation is blocked in 
Type I explosions by metal lines. Contamination from SN 
Ib/c can be particularly relevant at hi gh redshift and affect 
the cosmological parameter estimation (Homcier 2005). Effi- 
cient photometric identification of Ib/c is also possible using 
four-band photometry focused on the bluer bands a nd with 
modest photometric accuracy JCal-Yam et al. 20041) . At the 
highest redshifts (z > 1) weak lensing from the intervening 
mass distribution introduce a scatter of about 10% (with 
zero shift in mean), and even weak lensing maps can cor- 
rect only a small fraction of this effect (Dalai et al. 2003). 
This scatter is much smaller at lower redshifts (z < 1) and 
can eventually be inferred from the data (Amanullah, Mort- 
sell & Goobar 2003). Due to the large number of SN la per 
redshift bin this effect is negligible and averaged out. 

A sample of tens or hundreds of SN la studied in this 
way might reveal the spectroscopic characteristics associated 
with these photometric properties. Furthermore, if we can 
also measure the redshifts of these SNe's host galaxies, we 
can derive accurate luminosities. 

Although ALPACA is conceived as a drift-scanned 
imaging telescope, it is not impossible to suggest how a spec- 
trograph might be introduced at a later time. This instru- 
ment would also drift scan, this time through the reimaging 
system of the spectrograph onto the time-delay integration 
(hence drift-scanned) CCD detector at the back of the spec- 
trograph. A spectrographic slit mask could also be scanned 
across the focus of the telescope in front of the spectrograph 
in a way which might deliver many spectra simultaneously. 
The implementation of this idea might also deliver spectra 
for many thousands of supernovae, and could study huge 
numbers of galaxies as well. 

Measuring the apparent brightness and host galaxy red- 
shift for these supernovae would allow us to compute their 
luminosities to the level of 1% or better relative to other SN 
within redshift intervals of Az ~ 0.1 (even given the current 
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uncertainty in the background cosmology) . This will give us 
the opportunity to isolate samples of Type la SNe within 
redshift bins of Az ~ 0.1. In one year's survey, each bin will 
contain thousand of SNe. Comparing these large samples of 
SNe of nearly identical redshifts (therefore nearly equidis- 
tant), will allow us to seek for subclasses of SNe la within 
each redshift bin. Such subclasses could be distinguished on 
the basis of temporal evolution in various bands (for ex- 
ample different expansion rates at the same peak luminos- 
ity would indicate different temperature evolutions beyond 
those implied by 56 Ni mass), or from higher-order moments 
of the light curves. If the fractional contributions of the var- 
ious subclasses vary with redshift, we could then compen- 
sate for the effect. Furthermore, the various subclasses will 
be checked against one other to identify and correct for any 
systematic drifts within individual subclasses. 

3.2 Baryon Acoustic Oscillations 

Baryon oscillations in the galaxy power spectrum have 
recentl y emerged as a promis i ng probe of dark en - 
ergy ka,ke fc Olazehrook 200.1 ISeo fc Eisenstein 2003h. 
The large-scale linear clustering pattern contains a series of 
small-amplitude, roughly sinusoidal, modulations in power 
of identical physical origin to the acoustic peaks observed in 
the CMB. These oscillations encode a characteristic scale, 
the sound horizon at the drag epoch, which can be used as 
a cosmological standard ruler in the low-redshift Universe. 
When applied in the tangential and radial directions, this 
method enables the measurement of the angular diameter 
distance Da(z) and Hubble parameter H(z) in units of the 
sound horizon, over a series of redshift slices, which in turn 
may be used to place accurate constraints on dark energy 
models. 

The preferred scale was recently identified in the clus- 
tering patter n of Luminous Red Gal axies in the Sloan Digital 
Sky Survey l|Eisenstein et al. 20051) . constituting an impor- 
tant validation of the technique. A power spectrum analysis 
of the final 2-degree Field Galaxy Redshift Survey produced 
consistent results (Cole et al. 2005). The baryon oscillations 
method is robust against systematic errors. The underly- 
ing structure of acoustic peaks can be modeled very accu- 
rately using the linear physics of the CMB. In the galaxy dis- 
tribution there are modifying "non-linear" effects owing to 
redshift-space distortions, halo bias and non-linear growth of 
structure. However, these effects are confined to broad-band 
changes to the power spectrum and do not introduce pre- 
ferred scales that could be confused with the early- universe 
sound horizon. Recent simulations have confirmed that the 
baryon oscillation signature is preserved with only minor 
degradation (Seo & Eisenstein 2005, Angulo et al. 2005, 
Springel et al. 2005, White 2005). The acoustic signature 
may also be measured from photome tric re dshift surveys 
l|Seo fc Eisenstein 200.4 iBlake fc Bridle 2004h . although the 
smearing of radial information implies that only Da(z) can 
be determined with any confidence. However, given a suf- 
ficiently deep survey and accurate photometric redshifts, 
baryon oscillations may be detected with only ~ 1000 deg 2 
of imaging in broad redshift slice s centered at z = 1 and 
z = 3. We use the technique of ijBlake fc Bridle 20041) to 
simulate the resulting measurements of Da(z = 1) and 
Da{z = 3). The expected galaxy number guarantees that 



the clustering measurements are limited by cosmic variance 
and not by shot noise (nP = 3 in the notation of Seo & 
Eisenstein 2003). This is entirely consistent with the antic- 
ipated survey depth of ALPACA. We assume a survey area 
of 775 deg 2 . For some galaxy classes, accurate photometric 
redshifts are p ossible with five broad ban ds, for Luminous 
Red Galaxies iPadmanabhan et al. 20051) cr z = <jo(1 + z) 
with <7o = 0.03 at z < 1, due to the strong spectral break 
which implies rapidly c hanging colours with redsh ift. For 
other classes <to = 0.06 iFernandez-Soto et al. 20021) . while 
good photometric performance should also be possible for 
z — 3 galaxies using the Lyman break. For a recent analysis 
see liber et al. (2006). We suppose that these accuracies are 
achievable over the ranges 0.5 < z < 1.5 and 2.5 < z < 3.5. 
(We omit 1.5 < z < 2.5 given that photometric redshifts will 
not be accurate in this range in the absence of near-infra-red 
data). For comparison, we also display the equivalent mea- 
surements from a spectroscopic survey with the same con- 
figuration, using the method of Glazebrook & Blake (2005) . 
Our results are listed in Table [5] An additional potential 
systematic error for an ALPACA large-scale structure sur- 
vey is the survey geometry, which is a long, narrow stripe. 
Statistically, the number of measured Fourier modes is in- 
dependent of the survey area. However, in order to robustly 
measure the clustering amplitude on a given scale, that scale 
must significantly exceed the minimum dimension of the sur- 
vey. Otherwise we risk an enhanced cosmic variance owing 
to our consequent inability to measure the underlying large- 
wavelength modes in that dimension. In our case, the stripe 
width (2.5 deg, equivalent to a projected 101 Mpc/h at z = 1 
and 194 Mpc/h at z = 3) is comparable to some of the scales 
of interest for baryon oscillations (20 — 100 Mpc). Further 
simulations are required to quantify the level of risk for en- 
hanced cosmic variance. 

3.3 Weak lensing and other datasets 

ALPACA will also provide a very deep imaging survey per- 
fectly suited for gravitational lensing studies. After three 
years of the survey, the cumulative AB magnitude limit 
in the r-band will be ~ 27, providing > 40 gal/arcmin 2 . 
The imaging quality will be limited only by seeing at the 
CTIO site, thus no worse than 1 arcsec. With this resolu- 
tion ALPACA will be capable of mapping the dark matter 
distribution from cosmic shear measurements with unprece- 
dent resolution. Small scale resolution (on the arcmin scales) 
combined with measurements up to 1° are particularly im- 
portant to constrain dark energy. For a given amplitude of 
the scalar perturbations the balance between the non-linear 
scales an d the large ones is sensitive to the dark energy pa- 
rameters l|Benabed fc Van Waerbeke 20041 ). A more p ower- 
ful method is to use lensing "tomography" <Hu 1 009) from 
different redshift bins. ALPACA's sensitivity ensures that 
the noise coming from the intrinsic ellipticity of galaxies 
is negligible at all scales, offering an optimal measurement 
of gravitational lensing in the range of interest. A known 
source of limitation arises from the correction of galaxy 
shapes for Point Spread Function (PSF) effects. Foreground 
stars in the ALPACA stripes can be used to determine the 
anisotropic PSF and correct the images with standard tech- 
niques (Kaiser, Squire & Broadhurst 1995; Heymans et al. 
2005). In contrast the isotropic PSF tends to circularize the 
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galaxy shape and therefore to dilute the lensing signal. In 
such correction factor proportional to the seeing 

must be applied. A seeing < l" is necessary in order to 
keep the uncertainty o f the shape correction factor below 
5% dErben et al. 200 J) . Improving the observational accu- 
racy is not the only complication, on the theoretical side 
we have a limited knowledge of the growth of the non-linear 
structures. Using N-body simulations it has been possible to 
predict the non-linear power spectrum u p to 1 arcmin with 
5% accuracy at most JSmith et al. 20031) . 

An interesting characteristic of ALPACA is that the 
same region of the sky will be sampled ~ 1000 times. This 
is an entirely new regime of the PSF correction, which has 
never been explored before. In fact the number of dithered 
images with current lensing survey is ~ 10. The larger num- 
ber of exposures could greatly help to better determine the 
effect of systematics. We leave the analysis of ALPACA weak 
lensing observations to future studies. 

The survey will measure the location and collect pho- 
tometric redshifts of ~ 10 s galaxies over a large range of 
redshifts, potentially leading to the identifications of many 
galaxy clusters and groups. This will allow for an accu- 
rate number counting test, which in combination with exist- 
ing measurements of Cosmic Microwave Background power 
spectrum can provide com plementary co nstraints on the 
dark energy parameters iWang et al. 20041) . Another use of 
ALPACA would be to measure the cross-correlation between 
the survey and the CMB maps for measuring the Integrated 
Sachs-Wolfe (ISW) effect (Crittenden & Turok 1996). Al- 
though due to the limited sky coverage of the survey these 
measurements may not be able to pr ovide competitive con - 
straints on dark energy parameters iPoeosian et al. 2 005). 
The survey will also provide a large sample of quasars 
for Alcock-Paczynski test (Calvao, De Mello Neto & Waga 
2002), as well as strong-lens time delays (for 10-20 QSOs) 
JShapiro 1964) . 



4 FISHER MATRIX ANALYSIS 

We use the Fisher matrix approach to forecast constraints 
on dark energy parameters from ALPACA observations of 
SN la and acoustic baryon oscillations. The Fisher matrix 
provides a simple and practical method for estimating pa- 
rameter errors. It gives a local approximation to the likeli- 
hood surface about a fiducial model for a set of independent 
observations (Tegmark, Taylor & Heavens 1997). Formally 
the Fisher matrix reads as 



m(z) = M +5\ogV L (z), 



(2) 
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where 9^ are the model parameters, Oi are the observations 
of the quantity O and (j{Oi) the experimental errors. The 
derivatives in Eq. Q are computed at the fiducial model 
parameter values. The covariance matrix is given by the in- 
verse of Fisher matrix, C = F~ x , with the one sigma errors 
on the model parameters given by the diagonal components, 



a e . In case of multiple independent datasets the 



combined constraints can be inferred by simply adding the 
Fisher matrices of each experiment. 

Supernova la observations measure the luminosity dis- 
tance dh through the standard-candle relation, 



where M = M — 5 log Ho + 25 is the "Hubble-constant-free" 
absolute magnitude and T>l(z) = HodL(z) is the "Hubble- 
constant-free" luminosity distance. In flat universe the lu- 
minosity distance is given by 

d L (z) = (l + z) f^4, (3) 



o H(z')' 



with 



31 1/2 



H(z) = H [(1 - n m )f DE (z) + Q m (l + z) 
where Q m is the dark matter energy density and 

r f \ /q f Z 1 + w ( z ') J ' 

fDE(z) = exp t 3 / l + dz 



(4) 



(5) 



with w{z) being the dark energy equation of state. Sev- 
eral parameterizations have been extensively studied in the 
literature, for simplicity we li mit our analysis to dark en- 
ergy models parametrized by l|Chevallier fc Polarski 200 ll : 
iLinder 2003T) : 



w(z) = U)0 + W a 



1+Z 



(6) 



We have a three-parameter space consisting of fl m , wo 
and w a . In addition the magnitude-luminosity distance rela- 
tion depends on the off-set parameter M which we marginal- 
ize over. For the fiducial cosmology we adopt a model with 
fim = 0.3, wo = — 1 and w a = 0. We compute the SN la 
Fisher matrix using Eq. Q and the derivatives of di(z) 
with respect to the model parameters using double-finite 
differences about the fiducial model. 

We account for the uncertainty of the observed magni- 
tude in each redshift bin, using the standard approach of 
adding in quadrature the statistical and the irreducible sys- 
tematic errors (see for instance Kim et al. (2004)), 
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with o s tat = 0.15 and JV; being the number of SN la in the 
i-th bin. Assessing the amplitude of the expected level of 
systematic uncertainty a syB is more difficult. Because obser- 
vations at higher redshift are more challenging the system- 
atic error should increase, this increase is usually assumed 
to vary linearly with the redshift, 



Cays = Sm- 



Zma J 



(8) 



where 8m is the amplitude of the systematics at the maxi- 
mum redshift bin of the As we have discussed 
in Section r3.ll we still lack robust quantitative prediction for 
SN la light curves and spectra. Although results from nu- 
merical simulations remain controversial, they seem to sug- 
gest that nightly high signal-to-noise multiband SN la light 
curves may suffice to reduce systematic uncertainties to few 
% level (see for instance Hoflich, Wheeler & Thielemann 
1998). The use of multiple methods on a vast SN la sample 
such as the one provided by ALPACA is likely to reduce the 
level of systematics to 5 — 7% or better. Constraints from 
future SN la survey usually assume a more optimistic 2% 
uncertainty. We forecast constraints for both 5m = 0.05 and 
0.02. 

As discussed in Section ^, 21 baryon acoustic oscillations 
in the matter power spectrum provide a standard ruler which 
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Table 2. Expected uncertainties on x = (1 + z)cIa(z)/s and 
x' = H~ 1 (z)/s in different redshift bins, for different redshift 
measurements . 



z 


a x (cr = 0.03) 


a x (spect-z) 


cr x > (spect-z) 




1.0 


0.039 


0.013 


0.022 


0.5 


3.0 


0.020 


0.009 


0.016 



can be used to infer cosmological distance measurements. 
Oscillations in the tangential component of the matter power 
spectrum give the comoving angular diameter distance at 
the redshift slice of the survey in units of the sound horizon, 



x{z) = (l + z) 



d,A{z) 



(9) 



where cIa(z) = cIl(z)/(1 + z)' 2 is the angular diameter dis- 
tance and the sound horizon s is given by 

r^drag 

s - I — J". (10) 



I) 



\/Qm d + fir 



where adrag is the value of scale factor at the epoch of 
baryon-drag (baryons decoupling) and Q. r is the radiation 
energy density. The sound speed of the barotropic fluid is 



2 1 

c« = - - 



1 



3 i + !e« 



(11) 



Oscillations in the radial component give the Hubble rate in 
units of the sound horizon, 



x'(z) = 



(12) 



Because of the dependence on the sound horizon, the pa- 
rameter space is extended to fib and the Hubble constant 
h — fzo/100. Again we use Eq. (0 to compute the Fisher 
matrix components. The expected errors on x and x' from 
ALPACA measurements are listed in Table |2] 



5 RESULTS 

We now discuss the results of the Fisher matrix analysis. 
We list in Tableland Table [I] the la limits on fl m , wo 
amd w a from SN la with Sm = 0.05 and 0.02 systematic 
uncertainty respectively. We also quote the errors inferred 
from assuming the standard prior aa m = 0.03 and compare 
with those obtained combining baryon acoustic oscillations. 

Supernova data alone poorly constrain the dark energy 
parameters due to the degeneracy in fl m , for instance for 
Sm = 0.05, we find an m — 0.62, a WQ — 0.49 and w a un- 
bounded. These errors reduce by nearly a factor 2 for the 
case Sm = 0.02 (Table [IJ. On the contrary assuming the 
standard fl m prior we obtain a WQ = 0.2 and a Wa =1.1 for 
5m = 0.05, and a WQ — 0.11 and a Wa = 0.53 for Sm = 0.02. 

Constraints on dark energy from baryon acoustic os- 
cillations are also limited by degeneracies with other cos- 
mological parameters, in particular fi m and h. In addition 
there is a further uncertainty due to the dependence on fib 
which determines the size of the sound horizon at decou- 
pling. Assuming these parameters to be perfectly known we 
find a WQ — 0.17 and a Wa = 0.70 for spectroscopic redshift 
measurements. In figure we plot the la contours in the 
too — w a plane from SN la (solid line) and baryon oscilla- 
tions (dash-dot line). 



»« o 



-0.5 




-1.5' — 
-1.3 



-0.9 



-0.7 



Figure 1. 1-cr Fisher matrix contours on wq — w a . The solid 
line corresponds to SN la with 2% systematic uncertainty and 
3% Gaussian prior on f2 m , while the dash-dot line corresponds to 
BAO measurements with spectroscopic redshift accuracy assum- 
ing perfect knowledge of Q m and h. 



We are particularly interested in the SN and BAO com- 
bined constraints, so as to determine how the combination of 
the baryon acoustic oscillations can realistically improve the 
dark energy limits from SN la under minimal prior assump- 
tions. In such a case we only assume a Gaussian prior on the 
baryon density, an b = 0.006 consistently with WMAP/BBN 
results (Spergel et al. 2003; Hansen et al. 2002) and h to be 
fixed. 

In the case of photometric redshifts, the combination 
of SN la and BAO gives an m = 0.01, a WQ = 0.19 and 
a Wa — 0.96 for Sm = 0.05, and an m = 0.01, a W(] — 0.11 
and a Wa = 0.52 for Sm = 0.02. The use of spectroscopic 
redshifts leads to slightly better constraints, for instance we 



find an n 



0.003, 



= 0.14 and a Wa = 0.58 for Sm = 0.05, 
and an m = 0.002, a WQ = 0.09 and a Wa = 0.43 for Sm = 0.02. 
The level of detection of baryon acoustic oscillations as ex- 
pected for ALPACA will provide poor information on dark 
energy beyond that which comes from better constraining 
Q m to the level expected from future CMB experiments. In 
TableQ^we list the expected parameter errors from combin- 
ing SN and BAO with CMB limits from Planck, the inferred 
constraints are compatible with those obtained in (Linder 
2005). Including the shear weak lensing can certainly im- 
prove these limits, we leave this to a future study. 



6 CONCLUSIONS 

Uncovering the nature of dark energy is the new challenge of 
modern observational cosmology. Over the upcoming years 
lots of effort will be dedicated to improving luminosity dis- 
tance measurements of SN la standard candles and extend 
observations to complementary dark energy tests, such as 
baryon acoustic oscillations and cosmic shear weak lensing. 
In order to accomplish such a program several experiments 
have been proposed both from ground and space. 

Here we have presented the ALPACA survey and de- 
scribed the main datasets that the project is expected to 
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Table 3. l-<7 uncertainties on Q m , wo and w a from SN la assuming 5% systematic error with Gaussian Q 
prior and with different combinations of BAO detections. 





SN la (5m = 0.05) 


+a Qm = 0.03 


+D A (z) (<j =0.03) 


+D A (z) + H(z) (spect-z) 




0.62 




0.01 


0.003 


O~w 


0.49 


0.2 


0.19 


0.14 


'w a 




1.1 


0.96 


0.58 



Table 4. As in Table|3]with 2% SN la systematic error. 





SN la (8m = 0.02) 


+<rn m = 0.03 


+D A (z) (a =0.03) 


+D A (z) + H(z) (spcct-z) 




0.30 




0.01 


0.002 


O~w 


0.20 


0.11 


0.11 


0.09 






0.53 


0.52 


0.43 



provide. The bulk of the data is represented by high signal- 
to-noise, multiband, nightly SN la light curves. At the end 
of the three year run ALPACA is expected to observe 
> 100, 000 supernovae distributed in the range 0.2 < z < 1 
and several hundreds at lower redshifts. This huge dataset 
will be useful for finding possible subclass of SN la at differ- 
ent redshifts and correlations in the multiband light curves 
beyond the brighter-slower relation. This will allow for cali- 
bration of the standard-candle relation and reduction of the 
systematic uncertainties. 

In this paper we have forecasted limits on dark energy 
parameters from the expected redshift distribution of SN la 
and the detection of baryon acoustic oscillations in the mat- 
ter power spectrum. Although the accuracy of BAO mea- 
sureme nts cannot compete with ded icated surveys such a s 
KAOS (iDev fc Bovle 20031) or LSST jTvson fc Angel 200 ll) . 
combining them with SN la data will provide competitive 
constraints under minimal prior assumptions. ALPACA will 
also deliver several other datasets and will be particularly 
suitable for weak lensing measurements. 

ALPACA stakes out a regime in survey parameter space 
that is unique and useful, the intensive nightly, multiband 
monitoring of fields of order 1000 deg 2 . While in principle 
other surveys such as LSST, Pan-STARRS, JDEM or DES 
are instrumented for such intensive monitoring, in practice 
they are driven by other motivations such as all-sky tran- 
sient monitoring or weak-lensing mapping, devoting much 
less time than ALPACA to the 1000 deg 2 goal. Indeed, us- 
ing these like ALPACA would be inefficient, since ALPACA 
is extremely economical for this application compared to one 
of these more flexible instruments. Liquid-mirror techniques 
have now reached such a high level of quality that liquid- 
mirror telescopes may be used for astrophysics and cosmol- 
ogy. The ALPACA project will take advantage of this and 
provide high quality data at a fraction of the cost of standard 
8-meter telescopes or space missions. The vast possibilities 
of the ALPACA survey are still not fully explored and fur- 
ther investigation is needed. 
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